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1. Introduction 



Consider the usual Hamiltonian H = ti + Yli>j v ij-> describing a two-fragment col- 
lision a+b^c+d between iV nucleons, elastic or inelastic, with or without rearrangement. 
It is customary to define the prior (post) interaction V = Yliea jeb Vi h W = ^2iec jed Vi j)^ 
respectively, at the cost of considering as distinct the nucleons of the projectile a (resp. 
ejectile c) with respect to the target b (resp. residual target d) in the initial (resp. final) 
channel. Under this lack of antisymmetrization, one defines an initial wave x, which is a 
straight product of plane waves for a and b. Similarly one defines a non antisymmetrized 
final wave x' f° r c an d d. The transition amplitude is then introduced as the on-shell 
matrix element < x' I (V' + V'GV) \ x > 5 where G is the full N — nucleon Green function 
at the energy E of the collision. 

In order to take the Pauli principle into account, the correct approach, as described 
for instance by Austern[l], consists in i) retaining XiX' as non antisymmetrized, but 
ii) calculating the matrix element < x' I (VA + V'AGV) | x >, where the projector 
A = (AH) -1 ) V V, namely the parity weighted sum of all the permutations V of N 

nucleons, takes care of the antisymmetrization. Note that the validity of the above for- 
mula for the antisymmetrized transition amplitude rests upon the assumption that initial 
and final wave packets x, x' are normalized to unity and describe non-overlapping frag- 
ments a, b and c, d, respectively. Asymptotic density and flux are then unchanged under 
antisymmetrization. We shall come back to this point in Section 5. 

The calculation of Born terms, antisymmetrized or not, being straightforward, the 
present paper fo cusses on the following generic problem: Let z be any complex energy 
(eventually z — > E + iO). Let Xa(^i---^N a ) be any wave function for the projectile, with 
internal antisymmetrization for these N a nucleons. For instance, Xa ma y be a boosted 
shell model wave function, with or without configuration mixing, but it can actually dif- 
fer arbitrarily from an eigenstate of the projectile Hamiltonian H a . Let Xb(fN a +i---TN) 
be an arbitrary, antisymmetrized wave function for nucleus b, and consider the straight 
product x = XaXb- For the final state, consider similarly x'ci Xd an d x' = XcX'di with 
antisymmetrization internal to c and d, separately. We want to calculate the number 

V(z) =< x' I V'A(z - H^V | X > ■ (1-1) 

The main motivation for this problem is, obviously, that V represents the multistep 
transition amplitude (T — matrix without Born term) when x, x' are channel waves and z 
reaches the retarded on shell limit. But the validity of the arguments listed in this paper 
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extends to all complex values of z and very general choices for x, x' ■ (The only restrictive 
condition for this validity is the square integrability of V | x > an d V | >, so that V 
be well defined as a finite number. It is automatically satisfied for two-fragment channels 
and short range potentials.) 

A variational calculation of V is described in Section 2. Then we introduce in Section 

3 a time independent mean field (TIMF) approximation [2] for this calculation. Section 

4 is a discussion of symmetries which can simplify the variational TIMF equations. A 
numerical application of the theory to a-a scattering is described in Section 5, making 
use of self consistent symmetries introduced in Section 4. The final Section 6 contains our 
discussion and conclusion. Finally we introduce in an Appendix a new formalism for the 
representation of nucleons in different clusters as fermions with different pseudo spins. 

2. Equivalent, Antisymmetrized Variational Functionals 

We will show that T> is the stationary value of either of the following three Schwinger- 
like functionals of two independent trial functions 



T 


=< | AV 


1 X > + < x' 


1 v | 


V > - 


< ^' 1 


(* - 


H)\ 




(2.1a) 


T' 


=< | V | 


X > + < X 1 


V'A\ 


^ > - 


< ^' 


(*- 


H)\ 




(2.16) 


F = 


=< tf' | AV | 


X > + < X 1 


V'A 


1 * > - 


- < 




-H) 


|*>, 


(2.1c) 



which differ by slightly different insertions of A. The variations of give the following 
stationarity conditions, 



AV | x > ~(z - H) | * >= and < x I V - < | (z - H) = for (2.2a) 

V | x > ~(z - H) | * >= and < X ' I V'U- < | (z - fT) = for JF', (2.26) 

| x > ~(z ~ H) | * >= and < x I VU- < | (z - H) = for F. (2.2c) 

For z complex, the full resolvent G = (z — if) -1 is a bounded, uniquely defined operator. 
Its action upon square integrable vectors V | x > 5 I X > 5 < %' I V' , < x' I 
returns also square integrable vectors. Thus stationarity is reached for a unique pair of 
trial functions, 

| * >= GAV | x > and < |=< x' I V'G for (2.3a) 
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I $ >= GV | X > and < |=< X ' I ^'.AG for JF', (2.36) 
| ^ >= GAF | x > and < |=< X ' I ^'AG for F. (2.3c) 
The corresponding stationary value of each functional is then 

T = (< X ' | V"G).AV | X > + < x' I F'(G^ | x >) - (< x' I - tf)(GAF | x >), 

(2.4a) 

F = (< x' I | x > + < X' I ^(GF I x >) - (< x' I - H)(GV \ x >), 

(2.46) 

F = (<x'\ V'AG)AV | x > + < X' I ^'A(GAF | x >) 

-(< x' I - #)(GAV | x >), (2.4c) 

and the idempotence of A, together with its commutation with H and G, reduce all nine 
terms in the right-hand sides (r.h.s.) of Eqs.(2.4) to the same value, namely V, see Eq.(l.l). 
Hence T = F = F = V. 

We notice that six out of these nine terms are integrals confined to an interaction 
volume, defined by the ranges of V | X > 5 < x' I V'. For z reaching the on shell limit, 
^5^' are no more square integrable states, but rather purely outgoing or ingoing waves. 
The calculation of T>(z), however, by means of any one among these six terms, and the 
calculation of its on shell limit V(E + iO), can be confined to a truncation of ^' to this 
interaction volume. The remaining three terms contain both \1/ and and can not be 
used to calculate T>(z) in the on-shell limit. 

For this reason, we now investigate an approximation in which ^l/,^/' are replaced 
by products <E>, $' of independent, square integrable, single particle orbitals (pi,^, re- 
spectively, $ = rii^i' ®' = iLv 7 ^ or antisymmetrized products of such orbitals (Slater 
determinants). According to Eq.(2.3a), ^' contains only the same partial antisymmetriza- 
tion as that contained in x' ■> and its approximation should thus read as a product 
&' c $' d of separate Slater determinants for the ejectile and the residual target degrees of 
freedom, while \& is fully antisymmetrized, hence $ must be a Slater determinant. This 
approximation for T thus reads 

T =< $' c $' d | AV | x > + < x' I V | $ > - < & c & d | (z - H) | $ > . (2.5a) 

In the same way, the corresponding approximation for F reads 
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F =< & | V | x > + < x' I V'A | $a$b >-<&\(z-H) \ $ a $ 6 >, (2.56) 

with partial antisymmetrization restricted to a and b. Finally for F both ^' are fully 
antisymmetric, and the approximation for F reads 

F = C + C'-D, withC=< | V | x >, C =< x I V | $ >, D =< | (z-lf) | $ >, 

(2.5c) 

where the insertions of A present in Eq.(2.1c) are useless since are Slater determi- 

nants. 

We now notice that the first term in the r.h.s. of Eq.(2.5a) and the second term in the 
r.h.s. of Eq.(2.5b) contain an explicit A, which projects products of Slater determinants 
$ a $b, $c^d> i n t° f uu determinants respectively. The same property belongs to the 

third term in the r.h.s. of Eq.(2.5a), because the implicit A contained in $ commutes with 
H and carries full antisymmetrization further, into the bra. A similar argument holds 
again for the third term of the r.h.s. of Eq.(2.5b). It can be concluded that JF, T' ', F, which 
are equivalent when are fully flexible trial fuctions, are still equivalent when one 

attempts the separation of single particle degrees of freedom. 

The next Section, Sec. 3, thus refers only to the variation of F, as defined by Eq.(2.5c), 
with Slater determinants <£' as trial functions. For obvious technical reasons, we will 
investigate a slightly less general problem than that announced when introducing Eq.(l.l): 
only those cases are retained where x (resp. x') is the product of two Slater determinants 
Xa, Xb, (resp. x'ci x'd)- The wave functions Xai---x'd f° r ^ ne nuc l e i involved in the collision 
correspond to boosted shell models, like those which define the initial states of TDHF 
calculations [3]. This technical restriction still defines a generic problem, however, since 
any matrix element of the operator V'AGV between correlated internal states of these 
nuclei can be expanded in matrix elements between uncorrelated states, via usual, static 
configuration mixings. 

3. Time Independent Mean Field Approximation 

Let {xai ol = 1, ...N a } be the set of (boosted) orbitals which describes the projectile. 
Let {x/3, P = N a + l, ...N} , { X ' 7 , 7=1, -N c } , {x' s , S = N C +1, ...N} be similar sets for 
nuclei b, c, d. It is convenient here to define the following wave functions, 
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N a N N c N 

i x > = n i > n i >, < % i = n < ^7 1 n < i> ( 3 - ia ) 

a=l /3=AT a + l 7=1 5=7V C + 1 

which are pure products of orbitals, and 



| $ >= det{| v?i >}, < |= det{< ^ |}, (3.16) 

which are determinants of orbitals, without the usual (AH) -1 / 2 normalizations, because it 
is slightly easier to calculate the contributions C, C to the functional F as 



C = {N a \N h \/N\f/ 2 <&\V\x>, C = {N c \N d \/N\) 1 ' 2 <tf\V'\Q>. (3.1c) 

As usual, we consider the matrices of the overlaps of the single particle orbitals, namely 
the matrices < ip\ \ Xj >, < x\ I <Pj >? < <&[ I <Pj > an d their inverse matrices Aij, A\-, Bij, 
respectively. This generates the standard cofactors M^, Mij knrn i,...M' i -,...N i j knm i 
such as defined by 



< | x >= det (< v?- I Xj >) = J] < ^ I > wi th My =< $' | x > A,-*, (3.2a) 

i 

= ^ < X5 I <Pi > M' ijlk with M^ fc =< X' I $ > - A kj A' u ), (3.26) 

NijmZ = X < </4 | ipn > N ijknrn i with N ijknml =< $' I $ > X 

(BuB m jB nk — BiiB mk B n j — BijB mi B nk +BijB mk B n i+ Bi k B mi B n j — Bi k B m jB ni ). (3.2c) 
Accordingly, one obtains the derivatives 

d < I $ > 

^ , — = Na, (3.3a) 

5 9M l kl = M' (3.36) 

a < x'i I ^ > 3 

dMklnm j.* in o \ 

— — — — — = MkUjnm- {6.6C) 
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A straightforward calculation then gives 

<& \ V\x>= g J2< <Pi<Pj I v I > Afi^a, (3.4a) 

ijaf3 

<x'\V\$> = \ E < *W I u I <W > M^ Sji , (3.46) 
< | (z - #) | $ >= zdet(< v?- I >) - 



^2 < <Si \ t \ (pj > - - < v'iV'j I v I WW > Nijik, (3 .4c) 

where Latin indices run from 1 to AT, while Greek indices a, [3, 7, 5 run only inside the 

labels allowed for nuclei a,b,c,d, respectively. All matrix elements of v are antisymmetrized. 

The functional derivatives of C and C then read 

<)C = (Njmi/M) 1 / 2 x 



8<4 



\ jap jmnctfl 



'iiiu,ij(3a I Xj ^ j ) 

(3.5a) 



5 I <Pj > 



(NMl/Nl) 1 ' 2 



x 



V ^7(5 ijSmn J 

(3.56) 

This suggests the definition of the following non Hermitian, channel mean field potentials, 

<-\S b \ X a > = < I X > _1 J2 < • ^ I U I > ( 3 - 6tt ) 

J/3 

< . 1 s a I x/3 >=< I x > _1 J2 < • ^ I v I > ( 3 - 66 ) 

< X ; I I • >=< x' I * > _1 E < \v\.(pi> M' Si , (3.6c) 
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< x'sl S'c I ■ >=< X\$>- 1 J2 < XsXi \ v\-<Pi> (3-6d) 

which correspond to traces of v on "densities" related to nuclei b,a,d,c, respectively. It will 
be noticed that St, acts upon orbitals Xa of nucleus a only. In the same way, S a , a trace 
on a "density" related to nucleus a, acts upon orbitals xp °f nucleus b only. Accordingly, 
both formulae, Eqs.(3.6a,b), can be interpreted as the definition of a unique mean field 
potential S. Completely similar considerations are valid with Eqs.(3.6c,d), both traces on 
nuclei d and c defining a potential S', acting in turn on orbitals x'-y an d x'51 respectively. 

The functional derivatives of D are 

D — = Y^ z - 1) I (pj > Nij -J2< ( Pk\t\<Pl> N klj i \(fj>- 

3 jkl 



6«? i 

1 y 1 

9 Y < • I v I <PW > N iklj ~ J Y < VmVn ' V ' ^ k ^ 1 > N mnijlk \ ¥3 >, (3.5c) 



j'fcZ jmnkl 



52<<p'i\(z- t) N ij -J2 <( p'k\t\ipi> <(p'i\- 



ikl 



5 \ (pj > 

\ Y < V'tf'k \ v\-<Pl> N iklj ~ \ Y < V'mVn I V I <Pk<Pl > ^mmjtt < ^ I • (3.5<Z) 
ifcZ imnkl 

This suggests the definition of the following non Hermitian, mean field potential, 



< . 1 u 1 ^ >=< 1 $ > 1 < ■ & 1 v 1 <w > ^ ( 3 - 6e ) 

< ^ I U I . >=< &\$>- 1 J2< <Pi<Pk \v\.<pi> N kl . (3.6/) 



ki 

With such definitions of S, S\ U, the functional derivatives now read 



() °- = (N a \N b \/N\) 1 / 2 <^\x>x 



Y\ S + Y < ^\^\U > A A\Xj> A 3l - s £ j <v' k \S\xi> A u A jk I X j > 

j V ki J jki 



jkl 

(3.7a) 



= (NJNdl/Nl) 1 ' 2 < x' I $ > x 



<5 1 ¥>j > 



S' 



ikl 



SD 



J2lz-t-U-J2<Vk\ + \<Pi>Bi k ] 



J2<V>'k\(t + U)\ipi> B h B jk | ipj > 



jki 



(3.76) 



<Pj > Bji+ 



(3.7c) 



SD 



6 \ cpj > 



Y t B ji <( P , i \lz-t-U-^2<<p , k \(t+^)\(pi>Bi k \ + 

. i \ kl " ) 



ikl 



< ip' k \{t + U)\ipi> B u B jk < ip'i 



(3.7d) 



We now take advantage of the fact that a Slater determinant is invariant, except for 
an inessential change of phase and norm, under the linear rearrangement of its orbitals. 
Hence a linear rearrangement of orbitals <p>j can diagonalize one out the four matrices 

< <p>\ | (fj >, < (p^ | h | ipj >, < x'i I <Pj > 5 < Xi I I <Pj > • I n the same way, a linear 
rearrangement of orbitals (p\ can diagonalize one out of the four matrices < cp^ \ <pj >, 

< ip\ | h | <fj >, < <p[ | Xj >) < I ^ I Xj > • I n a representation where both overlap 
matrices < tp\ \ Xj > and < x'i I > are diagonal, these formulae reduce to 



SC 



6<<d 



<<p'k 


s 


Xk > 




Xk > 



x*>-X) 



<4 


5 


X l > 


<4 





Xj > 



(3.8a) 



5C 



8 \ ifj > 



= (NJNdl/M) 1 / 2 <x'\®>< x'i I > _1 x 



< X 



<x' fc 


S' 


<pk > 


<x' fc 


<fk > 



E 



<*5 


5' 


<pi > 


< Xi 


<fi > 



< Xi 



(3.86) 



where one finds that the reference energy for S, resp. 5", is < | V \ x > I < & I X > 5 
resp. < %' | V | $ > / < x' I $ >; see a l so Eq.(3.9). Alternately, in a representation 
where both matrices < ip\ \ <pj > and < cp^ | h \ cpj >, with h = t + U, are diagonal, the 
simplification reads 
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5D <$'!$> 



SD <$'!$> 



z 



<<p'k 






<^'k 





where one finds the reference energy for h, 



<v'k 






<^'k 


<Pk > 







$ > 




$ > 



and self energies for cp^ 



<¥>'k 


(* + 


f) 




<^'k 







1 >, 
(3.8c) 


1 <Pi > J 


< V'j 1 h 




< f'j \ ( fj> ) ' 




(3.8d) 


> 


(3.9) 



<4 


h 


<Pi > 


< V 


t 

i 


<Pi > 



m = z- E + ^< '1*^ (3.10) 



We notice that the unknown orbitals </?j are much less coupled to one another in 
Eq.(3.8c) than in Eq.(3.7c), and the same remark holds for the <£>'• when Eq.(3.7d) simpli- 
fies into Eq.(3.8d). Hence, among all the possible rearrangements which could technically 
simplify practical calculations, we choose those representations which simultaneously di- 
agonalize < <p' i \ <p>j > and < <p' i \ h \ <p>j > . The remaining four matrices will in general 
be non diagonal. Eqs.(3.8c), (3.8d) also have the advantage that they are close to the 
traditional, homogeneous Hartree-Fock scheme. 



The stationarity conditions, 



then read 



S(C-D) S(C'-D) 

J =0,-- ^ = (3.11) 



(NJNbl/Nl) 1 / 2 < $' 


X > < <Pi 


P>i > 




$ > 



(vi - h) i & > = v a o,/ v — x 



Yl ( s + J2 < v'k 1 2 1 xi > Aik ) 1 Xj > Aji ~ S < ^ i 5 1 ^ > i > 



(3.12a) 
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and 



< <p'j I (Vj -h) = 



(NcW/Nl) 1 / 2 <x' 


$ > < (p'j 






$ > 



X 



Y^A'^ < x 'i I (s' + J2<Xk I y I w >4fc) -Z) <Xfc I S' I ^ >^; fc < x\ 

. i \ kl " / ikl 

(3.126) 

A slight simplification of the r.h.s. of Eqs.(3.12) is further possible: because Xa and Xb are 
Slater determinants, a linear rearrangement of the orbitals Xa among themselves, and a 
separate rearrangement of the orbitals xp-> °an partly diagonalize the matrix < (p[ \ Xj >, 
but the simplicity of Eq.(3.8a) cannot be obtained under such a partial diagonalization. 
The same holds for separate rearrangements of x'-y and x'51 which can induce the simplified 
r.h.s., Eq.(3.8b), partly only. 

It will be interesting, in future numerical applications, to investigate whether the 
residual off-diagonal parts of < cp^ \ Xj > and < x\ \ <pj > can be neglected. If so, 
the corresponding approximation would create a representation where all four equations, 
Eqs.(3.8), would be almost compatible. This "one-orbital-partner" approximation would 
represent the many-body collision by single-particle transitions Xi ~ ¥ v'n Pi ~^ x\i with a 
conservation of the label i. 

For the sake of completeness, we recall here a result, shown earlier[4] in a slightly 
different context where both Xi x' are full Slater determinants rather than products of 
such determinants. The result is still valid for the present paper: 

Theorem: The solutions (pi, (p'j of Eqs.(3.12) are consistent with the ansatz that both 
matrices < (p[ \ <pj > and < (p\ \ h \ tpj > are diagonal. 

Proof: For i 7^ j, multiply Eq.(3.12a) by < ip'j \ and Eq.(3.12b) by | ipi >, respectively. 
Both resulting r.h.s. vanish, because the antisymmetries of and $ with respect to their 
respective orbitals induce 

. . 8C 5C , 

Variation of C and C with respect to < tp\ \ and | tpj > and subsequent scalar multipli- 
cation by < <fj I and | <pi > will for j 7^ i create two identical rows ( or columns ) in the 
determinantal structure of C and C . 

4. Self Consistent Symmetries 
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In practical calculations one often uses channel wave functions which have some sym- 
metry. Then one can choose trial functions having the same symmetry, without destroying 
the self consistency of the mean field equations, if some commutation relations are fulfilled. 
Using such self consistent symmetries could be called a "projection before variation" and 
may reduce the number of coupled mean field equations enormously. 

Let us first consider a unitary operator P acting in single-particle space, which relates 
the Xi's through 

P\Xr >= e i9i \Xa(i) > and < x-\P ] =< x'^e'^ , (4.1) 

where a indicates a permutation among the orbitals and 9i is an arbitrary i-dependent 
phase. Then one can choose 

P\ipi >= e^K w > and < ^|pt =< ^ (l) |e"^, (4.2) 

if the following commutation relations are valid 

[t, P] = , [v, P <g> P] = 0, (4.3a) 

and if furthermore the transformation P leaves the orbitals in the same fragment or trans- 
fers all orbitals from one fragment to the other, as we shall prove below. 

For the overlap matrices =< (p'i\(pj >, a.ij =< <p'i\xj > 5 a ij =< XilVj >j the 
unitarity of P implies 

Pij = e^-^PvdMj) , atj = eW-Wa^w , = e^ 9 ^ 9 ^ a' a{i)a{j) . (4.4) 

Next let us study the transformation properties of cofactors and inverse matrices. Writing 
the generalized cofactor expansion of det /?, which is invariant under Eq.(4.2), as 

S ki detp = J2PijN kj =Y / ^ iei ~ ei) Pa(i)a( j )N kj , (4.5) 

3 3 

and comparing with 

8a(k)a(i) det/? = ^2f3a(i)jN a ( k )j = ^2 Pv(i)v(j) N (r(k)a(j), (4-6) 
3 3 

one finds : 

Ni^e-W-WN^w, (4.7a) 
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and for the inverse of the overlap matrix (3 

B ij = e i ^-^B a{i)a(j) . (4.7b) 

In the same way one gets 

Ma = e-^-^M a{l)a{j) , M-j = e-^-^M^y 
Eq.(4.2) together with Eq.(4.3a) give 

< <Pi\t\<Pj >= e^~^ < Va^WaU) > ■ 

and 

< <Pi<Pk\v\<PW >= e l ^ +9l - di -^ < <p' a({) <p' a(k) \v\<p aU) <p a(l) > • 

From (4.8a,b) together with (4.7a,b) we find for the self consistent one-body 
and h 

< (Pi\U\<pj >= e l ^- 0i) < <p' a{i) \U\<p a{j) > and < ip'^ip, >= e l{B ^ B ^ < ip' a{l) \h\ip a{j) >, 

(4.9a) 

or extending to the full Hilbert space of single-particle motion 

[U, P]=0, [h, P] = 0. (4.96) 

Eqs.(4.9a) and (4.9b) together with Eq.(4.7b) will prove sufficient for the l.h.s. of the mean 
field equations (3.12) to be invariant. We will now turn to the inhomogeneous terms on 
their r.h.s.. 

To this end, let us define p (resp. p') as a matrix with elements equal to if the 
orbitals Xi an d Xj (resp. x\ an d x'j) are from the same fragment and 1 if they are from 
different fragments. Then if P satisfies 

Pij = Pa(i)a(j) andp^- = p' a{i)a{j ), (4.36) 

we have 

<• \[S,P]\ X i>=0, <X-\[S',P]\ ->=0, (4.10a) 

where 

S\ Xl >= (deta)" 1 ^^ < • <p' k \v\ X iXi > M kh (4.106) 

kl 

< x'i\S' = (deta') _1 E^ < X^xikH • W > M' M , (4.10c) 

kl 
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(4.7c) 

(4.8a) 

(4.86) 
operators U 



use having been made of Eq.(4.7c) and the analogue of Eq.(4.8b) for the channel potentials. 
To prove the self consistency of the ansatz (4.2), one has to show that 



5F 



5<<P'i\ 



5F 



(4.11) 



This is exemplified for Eq. (3.7c) where the symbol r\ is used for the z-independent part 
of the self energy 



Pdet/3 



^(z -h- r])B jt + ^2< <fk\h\(fi > BuBjk 

j jkl 



\<Pj > = 



e i9i det/3 



^2(Z -h- 7])B a{j)a{i) + ^2< (p' a{k) \h\(p a (l) > 5a(I)a(i)-Ba(j)a(fc) 
j jkl 



\<Pa(j) >, 



(4.12) 

using (4.7b) and (4.9a,b). This is the desired result because in the sum a(j) can be replaced 
by j for every summation index. The same can be done for the r.h.s. of the mean field 
equations (3.7a) and for Eqs.(3.7b,d), which completes the proof of the self consistency of 
the ansatz (4.2). 

The second case are antiunitary symmetries like time reversal. Let us recall the 
properties of antiunitary operators A, 



A(X\lp >) = \*A\ip >, 
<<f/\(A\<p>) = (< <f/\A)\<p>< 



and 



A^A = AA^ = 1 or < >= (< tp'\A^)(Atp >)] * . 
Let us consider a situation where 

A\xi >= e ldi W a{l) > and < X ^ =< X*{i)\e~ l9i ■ 
Then one can choose 

A\<fi >= e l0 *\ip' a{i) > and < =< <p u{i) \e- i6 % 

to fulfill the mean field equations provided A obeys the following relations 

At = t^A, (A <g> A)v = v^(A®A), 
14 



(4.13a) 
(4.136) 

(4.13c) 
(4.14) 

(4.15) 

(4.16a) 



and leaves the orbitals in the same fragment or transfers all orbitals from one fragment to 
the other such that 

Pij = PaUMi) ' P'ij = P*UM*)- (4.166) 

For Hermitian t and v, Eqs.(4.16a) reduce to commutator relations. From (4.13c) it can 
be seen that 

Pij = e^'-^^OVfi , ay = e l ^-^o4 0>(l) , a' l3 = e^-^a a{j)a{l) (4.17a) 

and 

deta = deta'. (4.176) 
Like in (4.5) and (4.6), one has for example 

S ki deta = Y^ aij M kj = S a(k)a(i) deta' = ^e^""^^^^^.)^, (4.18) 

3 3 

hence we can conclude 

My = e-^-^M' aU)a{l) , M-j = e-^- e ^M a{jMl) , (4.19a) 

similarly 

N i:i = e -*(*-'i)jV ff(j>(i) , B l3 = e^-^B aU)a(l) . (4.1%) 
For the matrix elements of t one calculates 

< ip'^ipj >= e'^-W < >, (4-20) 

while for v 

< VWMVW >=< <Pi<p' k \ \v{A®A)^{A®A)\<p m > = 

([< ip'M(A®A)^ v t [(A ® A)\cpj(pi >])* = e *(*i+**-«,-*«) < ^yj^djlul^w^cfc) > . 

(4.21) 

Therefore 

< >= e***-'') < V^I^M > , < V^l%, >= e*^"^) < ^^Ifcl^i) >, 

(4.22) 

whereas 5 is related to 5" 

< ^\S\xj >= eW*-™ < x'„ U) \S'\<P„ {i) > , < Xi\S'\ Vj >= e«*-W < <f/ a{j )\S\x*V > ■ 

(4.23) 
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Now we are able to express the commutation relations for S, S' and U, h with A as 

AU = U^A , Ah = h^A, (4.24a) 

< • \(AS\ Xi >) =< ■ \(S'U\ Xi >) , (< x 'i\AS')\ ■ >= (< X 'i\sU)\ • >, (4.246) 
and as in (4.12) we can conclude 



Adet(3 



^2(z -h- rj)Bji + E < fk\ h \fi > B u B jk 

j jkl 



\<Pj >-- 



I E(^* " ^ - + E < >* BfjB*^ \^>. (1.25) 

3 jkl 

Apart from a phase factor we recognize the l.h.s. of the complex conjugate mean field 
equation resulting from 6F/8\(p a ^ >. Doing the same for the r.h.s. proves the consistency 
of the ansatz (4.15). It is worth noting that the above proofs simplify if one uses the 
diagonal representation in channel-spin formalism (see Appendix). 

It is possible to generalize (4.1) or (4.14). We consider the case where the unitary 
operator P acts only inside the given set of orbitals Xi an d Xi such that 



>= E l j*Ixj > and < ^' pt = E < x'j\ L i 



(4.26) 



Then the following ansatz is self consistent 

p\<Pi >= E L ^i^' > and < ^' pt = E < 



(4.27) 



i.e. it fulfills 



8<^\ 5{<<p'i\P + ) 



provided the commutator conditions (4.3a) are satisfied together with 

[p,L] = [p',L]=0. 



(4.28) 



(4.29) 



If (4.28) is valid then P generates out of the i-th mean field equation the equation for the 
transformed orbital, which therefore can be left out of the actual calculation. The proof 
starts by noting that 

/% = E4^, (4.30a) 



ki 
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in short 

(3 = L^/3L, (4.306) 

which leads to 

[P, L] = , [N, L] = and [B, L] = 0. (4.31a) 
Here N is the transpose of the cofactor matrix N t j. By the same argument 

[M, L] = and [M', L] = 0. (4.316) 

Using (4.3a) one finds 

Udetp = Y,<-<Pk\v\-<Pi> N ki = Pi [Y,H L ln L ml<- V'M- <Pl>N kl jP, (4.32) 

kl \ kl mn / 

and applying (4.31a) results in 

U = pU (det/3)" 1 J2 N ^m <-<p , n \v\"Pm>jP = P ] UP. (4.33) 

\ mn J 

Thus one can list the following properties 

[U, P}=0, [S, P] = [S',P} = and [h, P] = 0, (4.34) 

which lead to 

[U, L] = , [S, L] = [§', L] = , and [h, L] = 0. (4.35) 

U, S, S' and h are defined as the matrices Uij =< tp'^\U\tpj > etc. With these relations 
the proof of (4.28) is simple. For the example (3.7c) the l.h.s. of (4.28) can be written as 



(det(3) ^[(z-h-fj) (LB) oi + (LBhB) 



j -i 



m >-- 



= (det/3) J2[(z-h-fj) (BL) J% + {BhBL) j% ] \^ >, (4.36) 

3 

where one recognizes the r.h.s. of (4.28). The same generalization can be done for anti- 
unitary transformations. On the other hand, the invariance of a Slater determinant under 
unitary transformations u of the orbitals, with detu = 1, makes it possible to choose 
orbitals u\i so that (4.26) transforms into (4.1). The conditions (4.29) tell us that L has 
to be a block matrix like 

<o ° 2 ) OT < 2 o)- < 4 ' 37 ' 
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while u is restricted to type a), because the Xi an d X^'s are only partially antisymmetrized. 
With L, li and l<i unitary there exists a unitary matrix which diagonalizes them. Choosing 
such a u brings (4.26) into the form (4.1), with the e l0i as eigenvalues. For case a) one finds 
a(i) = i while in case b) one has cr 2 (i) = i where a(i) and i refer to different fragments. 
In the next Section we use self consistent symmetries to reduce the 16 equations for the 
alpha-alpha system. 

5. The a-a system 

The elastic alpha-alpha collision has been studied by many authors [5] as a test bench 
of microscopic theories of collisions. As an application of antisymmetrized TIMF we con- 
sider elastic a — a scattering in forward direction. A central, spin- and isospin-independent 
interaction is chosen as a sum of a short-range repulsion and a middle-range attraction, 
both Gaussians 

v(r, ?) = V a exp(- |r ~ 2 r/|2 ) + V r exp(- ^ ~ ^ ) . (5.1) 

As potential depths V a , V r and ranges r a , r r we use the parameters proposed by Volkov [6] 
V a = -83.34MeV r a = 1.6fm V r = 144.86MeV r r = 0.82fm . 
The channel wavefunctions are described by boosted single nucleon Gaussians 

x (r, k) =7r -ip-i exp(ik-r- ^ ) , (5.2) 

where k is the single particle boost. Let fi,f5 denote the coordinates of the proton with 
spin up in the projectile and the target, respectively. In the same way, let tq, r*3, r 7 , r 4 , r 8 
denote the protons spin down, the neutrons spin up and the neutrons spin down, respec- 
tively. The 8-particle channel wavefunction is therefore given by 

X(fi, f 8 , k) = xiifi, k) x ... x x 4 (r 4 , k) x x 5 (f 5 , —k) x ... x x 8 (r 8 , -A;) . (5.3) 

Our problem has the following symmetries: 

a. Spin Isospin symmetries within each cluster imply 

Xi=X2 = X3 = X4 and X 5 = Xa = X7 = Xs- (5.4) 

b. Elastic scattering in forward direction, bra-ket symmetry: 

Xi(r) = Xi(r) for all i (5.5) 
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c. Forward scattering together with time reversal: 



Xi(* ! ) = X5(* s ), Xm = x'i(?) (5-6) 

d. Parity: 

Xi(r) = Xs(-r), Xs(r) = x'i(-r). (5.7) 

Following Section 4, we can choose trial functions which are consistent with our vari- 
ational equations 

V?i(r1 = V?5(-r), V ' b {r) = V ' x {-r) (5.8) 
^(f)=^ 5 (r), ^(r) = ¥>'i(r)- (5.9) 

Spin and isospin symmetries are assumed again, naturally. As described in Section 3, it is 
convenient to use a biorthogonal representation : If we rearrange the trial functions into 
parity eigenstates, such a representation can be introduced as 

<P+(f) = \J\ (<Pi(f) + <P5(r)) 
<P-(f) = \J\ (-Vi(r) + ^5(r)) 

^-(r1 = y|(Vi(r1+^ 5 (r1). (5.10) 

The symmetries and the spin-isospin independence of the chosen interaction reduce the 16 
coupled equations to only two coupled, non linear equations for <p + and (f>- which have to 
be solved self consistently. 

Because of our use of wave packets, zero point energies have to be taken into account. 
The zero point energy of a single nucleon is 

3 Ti 2 

<x'i\U\Xi>k=o= E = (5.11) 

where m is the nucleon mass. Note that E is also the center-of-mass zero point energy 
of an alpha particle. The binding energy of the alpha particle consists of the internal 
nucleon-nucleon interaction and three zero point energies, 

E a = 3£ + 6 <x'iX>\xiX2> ■ (5.12) 
19 



We adjust the width of the wave packet to 1.4 fm to get the experimental binding energy 
of the alpha particle. With wave packets the relative motion of the two alpha particles is 
not a pure plane wave, but also a wave packet of Gaussian shape. For that reason, a zero 
point correction for the relative motion is necessary. 

Because of the Galilean invariance of our many particle system, the unphysical kinetic 
energy of the center-of-mass motion has to be subtracted. Hence the relative motion 
Hamiltonian reads 

1 8 

H r = H - T CM = H - — - Vp! 2 . (5.13) 

16m 

i=i 

The problem of Tqm subtraction is the occurence of non diagonal terms , such as pi ■ p 2 . 
It is much simpler to approximate the center-of-mass kinetic energy operator Tcm by its 
zero point expectation value, 

<T CM >=E . (5.14) 

With all these arguments, the total asymptotic energy consists of the kinetic energy of the 
relative motion and the binding energy of the two alpha particles as well as the zero point 
correction of relative motion and the center-of-mass subtraction 

h 2 

E = (4k) 2 + 2E a + 2E . (5.15) 

4m 

The two coupled equations for (p + and (f>- are solved numerically, using a scaled Hamilto- 
nian and the channel wavefunction as starting points. With the help of the Hamiltonian 

#a = J> + A^ % - (5.16) 

i i,3 

we modify the stationarity conditions (2.2c) by introducing an additional term 

\AV\x > +( x ~ A Yl \vn\x>- (z-H x )\$>= 0. (5.17) 

i,jel,- ■ - ,4 
i,je5,-.-,8 

We start with a large scaling parameter A, when Eq.(5.17) implies | $ >~| x an d reduce 
this parameter to the physical value A = 1. We also start with complex energies E + iF 
and decrease V to zero to get the on-shell limit of <x'|AT|x>. A partial wave expansion 
for the mean field solution is useful to reduce the numerical expense. Choosing k along 
the z-axis, one can expand the channel wave functions in partial waves with just magnetic 
quantum number m = 0: 

rXi(0 = 5>i(r)*iV), (5.18) 
l 
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with 



,2 



3 3 A 1 7" 

Xi (r) =r7r-i / 9-'/(47r)'i , j I (Ar)exp(-— j) . (5.19) 
In analogy we consider the following expansions of (p+ and </5_, see Eqs.(5.10), 

r£+(r)= Vi+{r)Y?{r\ 

l even 

T<p-(f) = Y,<Pi-(r)Y?{r), (5.20) 

I odd 

where the summation can be restricted to even and odd numbers, respectively. E.g., the 
left hand side of our mean field Equation (3.12a) in partial wave expansion for (p + then 
reads 

n 2 d 2 h 2 i(i + iy 



v+ + 



2m dr 2 2mr 2 



<p l+ (r) - 3 U+ + {r)ipi>+{r) - 4 £ U„,~ (r)<p v+ (r) 
v v 

+ J2 U w + (r)W-(r), (5-21) 



where 



^>) = £^'(ooo) ^ r(r) ' (5 - 22) 



j j' 



JJ'L 




¥>j+( r ')W+( r ')A l ++» for (tti/) = (++); 
x { ^_(r')¥y-(r')/ra__, for (tti/) = (— ); (5.23) 

<Pj-(r , )<py+(r , )/ n —i for H = (-+)> 

C L (r, r') = V a J L (-4- ) exp(-^-) + K 7 L (^f) exp(-^-), (5.24) 

r a r a r r r r 

n_|__|_ =<</5' + |</5_|_ >, n =<</5'_|</5_ >, (5.25) 

and is the self energy corresponding to Note that in general is different from 
rj + because of the chosen representation. Here II denotes the modified spherical Bessel 
function of order L. Also note that in Eq.(5.21), when written with explicit parity coupling, 
only local potential terms appear. 

Corresponding to the expansion of the l.h.s. of Eq.(3.12a) an expansion of the source 
term is possible by using x instead of (p. A channel potential S^/ 1 acting on xi'i can be 
defined as 

^' 1 W = E^ / (ooo) 2 ^' 1(r) ' (5 - 26) 
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with 

S£\r) = f dr'C L (r, r>) £ £ jj'L ( J Q ^ V Vj+ (r') Xfl (r'). (5.27) 

The spin and isospin symmetry reduces the large amount of co factors of first, second 
and third order to only cofactors between nucleons with same spin and isospin. The con- 
vergence of our method is usually fast. On the way to the on-shell limit neither bifurcation 
nor crossing occurred. As described in [7], solutions with incoming waves sometimes creep 
in. Standard relaxation techniques [8] change these solutions into outgoing ones. Relax- 
ation techniques are also useful to speed up the convergence of our iterative method and 
remove other unphysical solutions. 

The Born terms of both the antisymmetrized (o) and the non-antisymmetrized (o) 
theories, respectively, are shown in Figure 1. Note their different behaviours when the 
boost k diminishes. This is an obvious illustration of the importance of exchange terms at 
low energies, whereas the Pauli principle may be neglected at high energies. One should, 
however, realize that the wave packets %, x' whose widths are related to the sizes of the 
fragments, may poorly correspond, at low energy, to well separated asymptotic fragments 
as required when using the above formula for the antisymmetrized transition amplitude. 
In consequence, the antisymmetrized wavepacket A\x > has vanishing norm as k tends to 
zero, and the antisymmetrization effect seen in < x'l^«4|x > is likely to be too strong at 
low /c-values: If the relative motion part of x were a plane wave, one would only have to 
exclude k = where the fragments fully overlap; since we deal with wave packets, we have 
to exercise care in a finite region around k = 0. This defect can be remedied completely 
by taking suitable superpositions [9] at different /c-values on both the bra and the ket 
sides in the transition amplitude. Rather than going into such an elaborate calculation we 
restrict ourselves to estimating the correction by dividing < x'|V.A|x > by the required* 
normalization factor a/ < x|-4|x > on the ket side. As a result the antisymmetrized Born 
term is modified at low k- values as shown by the (A) curve in Fig.l. The limiting value 
remains zero nevertheless. 

Figure 2 displays the real parts of <x'|AT|x> as a function of the total asymptotic 
energy E for both the antisymmetrized and the non antisymmetrized calculations. In the 
same way, Figure 3 shows the energy dependence of the relative difference A between the 
antisymmetrized and non antisymmetrized results for the imaginary part of the multistep 

* Note that V and A do not commute, < >^< x'\AVA\x >, and renormal- 

ization is not required on the bra side. 



22 



amplitude 

. < X '\lraAT\x>NA - <x'|ImAT|x>AS 



(5.28) 



<X / |ImAT|x>7VA 

As expected, antisymmetrization is important in the low-energy region. Regarding the 
defect of the channel waves x, x' at low energies, the same arguments apply as given above 
for the Born term. 

Finally we show on Fig. 4 an example of one of the wave functions cp generated by our 
code. It is remarkable that it is a square integrable state, although we are dealing with a 
theory of collisions. The reason for this is not only that the source terms x, x' are wave 
packets. Rather, a main feature of this formalism is the ocurrence of retarded self energies 
rj with finite imaginary parts even when the imaginary part of the many body energy E 
vanishes on shell. We notice on Fig. 4 the damped oscillations of away from the center 
of the interaction region. The damping is indeed due to the retardation expressed by Im 

V- 



6. Discussion and Conclusion 

At low energies not exceeding a few MeV per nucleon, previous microscopic theories of 
collisions [5] have been quite successful in the description of nuclear collisions. This success 
can be understood, to a large extent at least, from the limited number of channels open 
to the reaction mechanisms, as analyzed, e.g. by Tang [10]. Conversely, at large energies, 
beyond those corresponding to the Fermi momentum, effective theories have been also 
successfully proposed [11] for the derivation of optical potentials between the colliding 
nuclei. There is still room for a new microscopic approach, as shown by the present paper. 

Indeed we have proved here that the TIMF antisymmetrized theory can be imple- 
mented numerically, although it seems to be somewhat demanding and ambitious in ac- 
tivating all degrees of freedom. The TIMF equations involve all the single particle wave 
functions on the same footing, while we can insist that low energy [5] as well as high energy 
[11] theories essentially privilege the wave functions and/or the interactions governing the 
relative degrees of freedom between the clusters. These TIMF equations generate non per- 
turbatively intermediate states $, <&' which parametrize the multistep transition amplitude 
and give a picture of the reaction mechanism, as illustrated by Fig. 4. 

Once again we want to stress the fact that we are using a wave packet representation of 
the T-matrix rather than the more physical plane wave representation of relative motion. 
The conversion of one representation into another [9] is a typical generator coordinate [12] 
problem, which we intend to study. 
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Finally, it is clear that the mean field calculation which has been presented here can 
accomodate higher order corrections. An advantage of the TIMF technique is its close 
connection with the traditional techniques of the shell model. Therefore Tamm-Dancoff 
corrections and even RPA corrections are available. 

Some parts of this paper are also part of the Diplomarbeit of A. Wierling and the 
Doktorarbeit of J. Lemm at the University of Miinster, respectively. This work has bene- 
fited from a Procope grant for franco-german scientific cooperation. Two of us (B.G. and 
A.W.) thank the nuclear theory group of Kyoto University, where part of this work was 
done, for hospitality. 
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The prior and post interactions V, V are symmetric under nucleon exchange within 
the respective fragments but unsymmetric when nucleon exchange between fragments is 
included. Under this partial lack of permutation symmetry, only pure product or in- 
trafragment antisymmetrized wave functions x, x' were used so far in this paper. This is 
similar to the second quantized formulation of TDMF by Reinhardt [13]. We shall now 
show that there exists an equivalent "channel spin" formulation which works with totally 
antisymmetrized, two-component channel wave functions. For the sake of clarity of pre- 
sentation, we shall slightly modify our notation of Sections 2 and 3 in the following way: 
We denote the pure products of (normalized) single-particle wave functions in initial and 
final channels, eq(3.1a), by round brackets |x) and (x'|, respectively. Curly brackets are 
used for channel wave functions 

|X} = VNJN b \A a A b \ X ) , (Ala) 
{ X ^ = ^NjNJ.(x\A c A d , (Alb) 

which are normalized and antisymmetrized within the fragments a, 6, c, d, in obvious 
notation. Angular brackets, like in |<& >, are reserved for fully antisymmetrized wave 
functions, 

Ix>=^[EH^Ix)- (Ale) 

With such |x}, {x'l it is in general not possible to diagonalize all three overlap matrices 
a,ij =< (p^Xj > 5 <x'ij =< XiWj >i Pij = < f'iWj > by rearrangement of orbitals within 
the four sets of wave functions {xi}, {x'i}, Wi}, Wi}- This would correspond to solve the 
following set of equations for the matrices M , M' and B, B' where the latter two have 
block structures which do not mix orbitals of different fragments, 

M'I3M = dp , 

M'aB = d a , (A2) 
B'a'M = d a > , 

with diagonal matrices d a , d a >, dp. This means that B~ l a~ l f3a'~ 1 B'~ 1 = d^dpd^} has 
to be diagonal, which is not possible in general because of the block structure of B and 
B' , which is a consequence of partial antisymmetrization. To diagonalize all three overlap 
matrices one has to work with totally antisymmetrized channel wave functions x and x' ■ 
Because the actions of V and (H — E) on an exact channel eigenstate are equivalent, 
one could replace V by (H — E) which commutes with the total antisymmetrizer A. How- 
ever, the equivalence of (H — E) and V (or V) is destroyed if approximate channel wave 
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functions Xi x' are used, like in the mean field approximation, when {H — E)\x > and 
V\x > may differ by long-ranged terms. Because V is normally short-ranged, V should be 
used. In order to alleviate the clumsy partial antisymmetrization imposed by V and V', 
the concept of channel (pseudo-)spin is introduced in this Appendix. 

Inside the matrix elements which involve x or X* ■> the single-particle wave functions 
Xi, x'i ari d <fi, are translated into two-component (channel spin) wave functions 



\Xi >- 



\xT >-- 



\xt >^ \xT >-- 



<x'i\^< xr\ 



<x'i\^< x'r\ 



Xi 




Xi 



x\ 





Xi 



for iea, 

for ieb, 
for i e c, 
for ied, 



\<Pi >^ \<p? >-- 



<Pi 



Therefore we consider the transformation 



Xi 




jeb 




(A.3a) 

(A.Sb) 
(A.3c) 
{A.Zd) 
{A.Se) 
(A.3f) 



(AAa) 



{X\ - {X 



ICS I 



X'i 




\ 



jed 



(AAb) 



For x a channel spin 'up' means that the orbital belongs to fragment a, and channel spin 
'down' means that it belongs to fragment b. For x' the distinction is between fragments c 
and d. Equal components in channel spin formalism, eqs. (A.3e) and (A.3f), are assigned 
to single-particle orbitals (p i: (p[ of trial functions $, This assignment has no effect on 
total antisymmetrization, 



|$ >^ \$c* >= V \VU J| 



Pi 



{AAc) 



A. 



{AAd) 
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Now we introduce channel spin dependent potentials which replace V and V but leave 
their matrix elements invariant. With the Pauli matrix a l z acting on particle % in %-channel 
spin representation and with a z acting on particle % in x'-channel spin representation, the 
replacement of the interaction Vij between particle % and j by 

vfj = -(1 — o~ l z (j J z )vij in x-channel spin representation, (A. 5a) 

v'^ s = ^(1 — a' z l a'J)vij in x'-channel spin representation, (A.5b) 

sets the interactions vf?, v'^ s to zero if acting within fragments and leaves them equal to 
otherwise. With the definition 



tea 
jeb 



2 



zee 
jed 



2 



it is easily seen that 



< &\V\x} =< $ /cs |F cs |x cs }, (A6a) 

W\V'\Q >= { X ' CS \V' CS \$ CS > . (AM) 
In this enlarged channel spin space V cs and V' cs commute with A, 

[V cs ,A] = 0=[V' cs ,A], (A.7) 

and therefore 

< $' cs \v cs \x cs } =< & cs \AV cs \x cs } = 



y/NjNj. < & cs \V cs AA a A b \x cs ) = \J < & cs \V cs \x cs >, (A.8a) 
and correspondingly with (A. lb) 



/ TV 'TV J 

{ X ' CS \V' CS \$ CS >= \j c — { < x' cs \V' cs \$ ca > . (AM) 

We recall here that AA a = AAb = A. It is thus possible to work with fully antisym- 
metrized wave functions only. The functional F of Eq.(2.5c) then reads in the channel 
spin formulation 



r = < & cs \v cs \ x cs > +\/^§r L < x cs \v' cs \$ cs >-< &\(z-h)\$ >, 

(A.9) 
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where the channel spin is not needed for the < $'\(z — H)\G> > matrix element. (An 
introduction of channel spin in this matrix element is nevertheless possible and would 
require a normalization factor for every orbital cpi, cp^, but has no further effect. Note 
that this normalization factor must be included inside < $'\(z — H)\$> > only.) 

Within x cs and x' cs the orbitals can now be rearranged like those within $ and <&'. 
This leads to generalized channel spin orbitals 



\xT >= 



^ ) , (A.lOa) 



.. test _ I Xi 
^ Xi \ — \ 12 
Xi 



(A106) 



It is now possible to diagonalize the generalized overlap matrices < <^ cs |Xj S >; < x'^l^T > 
and corresponding generalized co factors M^, M[ c - S through linear rearrangement of the 
Xi & x'i S leaving the <fi, (p[ unchanged. This enables us to diagonalize < y^l^j > and 
< ip'^hlifj > by linear rearrangement of the cpi, (p[. In Eqs.(3.4a,b) the sum | Ylap nas 
to be replaced by \ Ylu where the indices k, I run over all orbitals because of the totally 
antisymmetrized wave functions 

< $'cs lV cs lx cs >= 1 J- < tf^T \v cs \x c k s xt s > M^ k , (A.lla) 

ijkl 

< x 'cs\ V 'cs^cs > = 1 J- < ^cs^cs^cs > M ,cs^ {A ub) 

ijkl 

Here antisymmetrized matrix elements are used for v cs and v' cs . 

The derivatives with respect to the single-particle orbitals (pi, ip\ must take into ac- 
count the fact that <fi resp. (p\ appear twice inside the channel spin trial functions. With 
the relations for generalized channel spin wave functions 

< vCvCK s |xrxr >=< f'mfMxlxl + Xlx\ >, (Al2a) 

< x'^XnWM 3 > = < x' m Xn+X m Xn\v\V>kW > (A126) 

for the antisymmetrized matrix elements, Eqs.(3.5a,b) become: 

6C =\[^^(lT,<^>\xlxl + xlxl>M^ k (AASa) 

V jki 

+ J E < tfmvMxlxi + xlx] > M™ nijlk \ x ) + x) > I , 

jmnkl / 
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Y, < ^? + x' k 2 x' l 1 \v\-^>M' k T jl 



(AA3b) 



jkl 



+ 7 E <xU? + XkX>\^n>M' k f jinm < X f + xf\ 



jmnkl 

Now the self consistent channel potentials can be defined by 

<-\Si\xl> 



< ^c! lx cs > E < • ^>\x\xl > M c n %, 



nk 



<-\s 2 \ x f> = 
<x' l 1 \s[\-> = 
<x?\s' 2 \-> = 



< &d\ x cs > E < ■ vMxtxl > M^ k , 



nk 



E < Xl l XkM ■ V?n > M^, 



< x' cs \& cs > 



kn 



< Y ,cLcs Z E < X?X*M • V* > M kn- 
X 1 kn 



(AAAa) 
(AA4b) 
{AAAc) 
{AAAd) 



Then Eqs.(3.5c,d) and (3.6e,f) do not change because there are no %, x' involved, but 
Eqs.(3.7a,b) now read 



SC = NJN b \ 
6<<d\ V Nl 



< $ /cs | x cs > 



x;^(5iix}>+5 2 i^>) 



+ ^ E (< ^ml-SlIX/ > + < Vml^lx? >) ^Ti A Tm\x) + x) > 



jml 



- E (< ^'m\Si\xl > + < ^'JS 2 \ X 2 >) A? m A%\ X ) + x] > 

jml 



, (Al5a) 



5C NANA 
6\<pi > ~ V AH 



E^S s (<x?i^+<xfi^) 



+ \ E (< xMbm > + < Xfl^bm >) < X? +X?\ 



jml 



~ E (< ^l^llVm > + < X? >) A'-^r < xj 1 + X 
jml 



(AA5b) 



These equations are invariant under linear rearrangement of the orbitals xf \ Xi \ an d 
they reduce to Eqs.(3.7a,b) when the wave functions (A.3a,b,c,d) are used. In a represen- 
tation where both overlap matrices < (p'f s \Xj S > and < x'^l^T > are diagonal, equations 
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(A.15a,b) reduce to 



SC 



NJN b l < $ /cs | x cs > 
N } - <ri\x}+Xi> 

1 



{s 1 \xl > +s 2 \xl >) 




< rtis^xl > + < ri\s 2 \ x f > \ , i 2 

/111 2 I \Xi ' Xi 

<<Pi\Xi+Xi> J 



-E 



<4 


Si\x] > + <( Pj 


S2\xl > 


<4 


x) + xj > 



(A16a) 
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x? + X? 


<pi > 
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<x? 


s'M> + < x ? 


s> 2 




< X? + X? 


<Pj > 



J2\ 



{A.m) 



Note that 



>= I! < M + X* > and < x' cs \® cs >= I] < X? + > • (^17) 



These equations can now be combined with Eqs.(3.8c,d) where < <p'A\pj > and < (p\\h\(pj > 
are diagonal. Because all the cofactors are simultaneously diagonalized when the three 
overlap matrices are chosen diagonal, the enormous number of cofactors, like iV 6 for 
M mn ijik with six indices, is reduced to the number of diagonal elements only. The price to 
be paid for this simplification is that one has to deal with two component wave functions. 
They give rise to additional terms in the inhomogenous parts of (A.16a,b) as compared 
to (3.7a,b). A generalization of this two-component channel spin formalism is possible 
for more than two fragments in one or both channels. For n fragments in a channel, n- 
component channel spin wave functions can be introduced in the corresponding matrix 
elements and give again a possibility of using full antisymmetrization. 
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